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The problem formulated concerns the mutual impedance between two vertical dipole 
antennas A and B located near the surface of a spherical smooth earth. The path between 
A and B is made-up of a number of homogeneous segments where the surface imj)edance is 
constant. Various formulas are developed, for two- and three-section paths, which are 
suitable for computation. Certain limiting cases are discussed and where possible a physical 
interpretation of the results is given. Comparisons with previous work are made. 

1. Introduction 

In tlie computation of the ground-wave field of a radio transmitter it is often assumed 
that the earth is a homogeneous and smooth conducting sphere. Wiiile this is a satisfactory 
approach for many applications, such as estimating the coverage of broadcast transmissions, 
there is often a need for more precise determinations. For example, the phase characteristics 
of the ground wave are significantly influenced by inhomogeneities in conductivity that occur 
at coast lines. This is an hnportant consideration in the prediction of errors in radio naviga- 
tional systems. 

An excellent semiempirical approach to the problem of calcidating fields over an inliomo- 
geneous earth has been given by Millington [1949]. Sometime earlier Feinberg [1946] had 
formulated a general method for treating mixed patlis over a flat earth. Similar results were 
obtained independently by Clemmow [1953] and I^remmer [1954]. The author [Wait, 1956a] 
showed that the integral formula for a two-section path could be easily derived by an appli- 
cation of the compensation theorem [Monteath, 1951]. It was also indicated in this paper 
[Wait, 1956a] that the results could be readily generalized to a spherical earth. In a further 
paper [Wait and Householder, 1957] extensive numerical results were given for propagation 
over a two-section path on a spherical earth. Furutsu [1956] has also considered propagation 
over a spherical earth under mixed-path conditions. His general method is based on solving 
the dual integral equations for the problem by an iterative procedure. 

It is the purpose of this report to discuss the theory for mixed-path propagation over 
both two- and three-section paths on a smooth spherical earth. Some attention is given the 
various representations which may be used in practical field computations. The situation 
when one terminal is at great heights is considered in some detail since this solution describes 
the radiation pattern of an antenna when the ground plane is inliomogeneous. Where possible, 
a physical interpretation of the results is given. 

2. Formulation for a Two-Section Path 

The mutual impedance Zab between two vertical antennas located at A and B over a 
spherical earth of radius a is considered. The situation is illustrated in figure 1 where a 




Figure 1. Two-section "path on a spherical earth. 



1 Contribution from the Central Kadio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. 
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vertical cross section of the earth is shown. The great circle distance between A and B 
(measured along the surface of the earth) is d. The earth medium to the left of the boundary 
line has a conductivity a and dielectric constant e. The corresponding constants for the medium 
to the right of the boundary are a^ and ei. A variable P on the great circle path is a distance 
a from B. Then, for a^d^, the surface impedance is 

Z^[^MoCo/(cr+^6a,)]^'^^l-^^T'' (1) 

while for a<Cdi 

L crii-2eicoJ 

where /xq is the permeability of the whole space which is assumed to be constant. 

The mutual impedance between the dipoles A and B for the two-section mixed path illus- 
trated in figure 1 is denoted Z\^. It was shown previously [Wait, 1956a] that it could be 
expressed in the form 

/2(Z;,-Z„,)^(Zi-Z) jj{H,,Hi,)dS (3) 

S 

where Zab is the mutual impedance if the surface of the earth were homogeneous with surface 
impedance Z everywhere. Hat is the tangential magnetic field of dipole A over the homo- 
geneous earth while Hbt is the tangential magnetic field of dipole B over the inhomogeneous 
earth. The current in the dipoles are both taken equal to / for convenience. The surface 
of integration S extends over the region of the earth which is characterized by a surface im- 
pedance Zi. Equation (3) also follows directly from the work of Monteath [1951]. 

3. Reduction to a One-Dimensional Problem 

Equation (3) is essentially a two-dimensional integral equation for the fields. On a 
straightforward application of the principle of stationary phase, the surface integral can be 
reduced to a line integral from a=0 to rfj along the great circle path. Before stating this result 
it is convenient to introduce certain attenuation functions, W and Wi, as follows 

Za,= ^4r^e-'''W(d,Z) (4) 



2Trd 
and 

2Td 



^ab 



e-'''W'{d,Z,Zi) (5) 



where la and /& are the effective lengths of the dipoles A and B. It is to be understood that 
W and W^ are also functions of ha and hi,, the height of dipole A and B above the ground. The 
functions W and W^ are normalized such that they would approach 1 if the earth were fiat 
and perfectly conducting — ^and provided Aa=A&=0. Consistent!}^ in what follows it is assumed 
that ha and hb<^d<^a. 

The resulting one-dimensional integral equation is given by 

W'(d,Z,Z,)^W{d,Z) forrfi<0 

l^ruJ'7 7^ w/^^N t^kd\'^f Z,-Z \ C'^ W(d-a,Z)W'(a,Z,,Z) , . , . ^ 
W\d,Z,Z.)^W{d,Z)^[^-^) (-^;J, [a{d-a)Yi^ rf.for^.>0, 

where r^o- W€o)^^'^120^. (6) 

This result was obtained in an earlier paper [Wait and Householder, 1957] where certain 
numerical results were given for a two-section path in the frequency range 20 to 200 kc/s. 
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In the present work, it is desirable to introduce certain diiiiensionless parameters. These 
are defined by 

iq: = ikal2y''(Z,M, 
x={kal2y"{d/a), 

x = (fa/2) 1/3 (a/a), 
Xi = {kal2y'Hdila). 
ThusW'(x,q,qi)^W(x,q)iovxi<^0, and 



W'{x,g,q,)^W(x,q)+l-j\ {qi-q)\ rV ^m./ ^■^ 

V"/ Jo \x(x—x)\ 



(V) 



for Xi>0. It is immediately observed that W^(x,q],q) in the integrand may be replaced by 
W(x,qi) since x'^Xi over the range of integration. Therefore 

W'ix,,,,o^W(x,a)+(^Y\a.-<i) p^^f^^l^^x. (8) 

\^^/ Jo [x{x—x)\ ' 

An alternate form of this equation is easily obtained b}^ regarding the left-hand portion of the 
path as a modification to the homogeneous earth of surface impedance Zi. Thus 

\^V Jo \x\x—x)\ 

4. Alternate Representations 

Equations (8) and (9) are explicit integral formulas to permit the computation of the 
attenuation function W in terms of the attenuation functions W appropriate for a homogeneous 
earth. The latter are well-known from the theory of van der Pol and Bremmer [1937] and 
Fock [1945]. From their work it follows that 

where 

x={kal2y'^ (dla), ya=(2/kay^'kha, y,= (2/kay^' kh,. 

The coefficients ts are solutions of the equation 

w[(t) — qwi(t) = (11) 

where Wi(t) is an Airy integral and the prime indicates a derivative with respect to t. In terms 
of Hankel functions of order one third, wM=exp (-2wi/3)(-7rt/3y^'Hi%{2/3){-ty/']. The 
notation used here is closely akin to that used by Fock [1945]. 
In a similar fashion, 



/I) 



where <"' are solutions of 

w[{t)-q,w,{t)=0. (13) 

Using (10) and (12), the integration in (8) may be readily carried out to yield 
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(14) 

The double summation converges quite rapidly if Xi is of the order of unity or greater and 
provided that either ya and y^ are not large compared with unity. 

When Xi is small, the double series expansion given in (14) becomes very poorly convergent. 
An alternate expansion for the case ht,=0 is obtained by using the following representation 
[Bremmer, 1953 and 1958; Wait, 1956]. 



m = 0, 1, 2, . . . 



where 



Using this result along with the residue series representation for W{x—x,q) enables (8) to be 
expressed in the form 

W^(x,q,q,)=W{x,q)-ix'^'{q^-q) ± ± A^e'^^'^q^Irn^ ^^lit.-y) ^^^^ 

where 

7^^= ^-z(.-.).,(^) 2 dx. (17) 

2 Jo 

The integral Im-i may be reduced by the following recurrence relation 

2 

4=^L-^(-^^)^s(5)-/2-.('-|^V^_1. (18) 

Successive application of this equation enables /(;;j_i)/2 to be expressed in terms of /q and /_i/2. 
These are given by 



(19) 



and 



'^fs 

/-i/2=6-^^^J^ -|— ^^=(^j 6-^-serf(V-^^.xO (20) 



where erf (Z) is the error integral of argument Z. 

The double series expansion given by (16) is highly convergent if x is somewhat greater 
than 1 and X] is somewhat less than unity. In fact, if |g?2i|<Cl, only the m=0 terms need be 
retained. Furthermore, if in addition x^<^\j 

And thus, 

W^{x,q,q,)^W{x,q) [l+^2 2^^} (21) 

The term in square brackets can be regarded as the correction to the attenuation function as a 
result of the portion of the path from J=0 to iCi. In terms of the original parameters of the 
problem, the square bracket term becomes 
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[l-^Y'' ^^-^ (2Wi)^/4 (22) 

This same correction emerged from the corresponding theory for the flat-earth case [Wait, 
1956b]. . 

5. Solution for One Antenna at Great Heights 

When one terminal becomes elevated to a large height such that ?/„ or y^^l, it is desirable 
[following Fock, 1945] to replace the attenuation function W by an auxihary function V. This 
function V actually characterizes the far zone radiation pattern of a vertical antenna on the 
curved surface [Wait and Conda, 1958]. First it is noted that if \y—t\'^l 

wi{t-y)^e-'^^\y-ty^ exp [-i ? (^/-O'^'T (23) 

Furthermore, if also y'^t, 

w,{t-y)^e-'''''y-'^'e-'h"\'y"'K (24) 

This latter asymptotic relation enables W{x,q) to be written in the form 



.2 ,/, /y.1/2 

W{x,q) ^e-' I "" Vix-^ly,.,i) ^ (25) 

ua 



where 



"(^•''--'"ScTi^^SytSr- « 



This suggests that the radiation pattern function V^(X,q,qi) for the mixed path be defined by 
the relation 

W'(x,q,qO = e-'l''''V'(x-^IV,,q,qd^,- (27) 

Ua 

The integral formula for V^ may then be obtained directly from (8), thus 

V'iX,q,q.) = nX,q) + ^f^ rYi^tzA^mmiJt (28) 

Kjf^) Jo (XY^^ 



where 



X=x-^'ya^y~) 



ka\^H--^2aha 



a 



It should be noted that X can be positive or negative depending on whether d is greater or less 
than {2ahay'\ 

An alternate form for V is obtained from (9) , thus 

where Xi=Xi — ^ya. Since ^a^l, the lower limit of this integral is effectively — oo . 

Series formulas for V'{X^q,qi) can be readily obtained in the manner discussed above 
for W\x,q,qi). However, these are applicable only in the present situation if X>0 since the 
series given in (26) is divergent for X<CO. A more suitable approach is to use the contour 
integral representation for V{X,q). This is given by 
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F(X-^)-2^^^ 



{t) — qwi{t) Wi{t) 



dt. 



(30) 



The contour is chosen so that it encloses (in a clockwise sense) the poles of the integrand at 
t=ts. It is readily verified that —27ri times the sum of the residues leads back to (26). It is 
now convenient to write 



ViXS^I^aMX" 



where 



^-inirl2 r 



r 



wi{t—y^) 
{t) — qwi{t)] Wi(t) 



dt. 



(31) 
(32) 



Methods for evaluating the coefficients a^(^) have been discussed by Logan [1959]. 

The above representation for V and the series form for the W given by (15) enable the 
integration in (28) for V^ to be readily carried out. The result is 



where 



Pn !^= f' {X-i)"ii) ^ dx. 
' 2 Jo 



(33) 
(34) 



The integral P^, !!izLi can be reduced by using the following formula 



(x-x,r(x^y 



n 



f+' 






(35) 



Successive applications enable P^ ^^ to be expressed in terms of P ^^ where 



2 Jo 



X - m — 



r''dx= 



{x,y^ 



2^ 



(36) 



The above double-series representation for V converges very rapidly when both \X\ and x^ are 
reasonably small compared to 1. 

6. Three-Section Path 

The next most obvious extension is to a three-part medium. The situation is illustrated 
in figure 2. The approach used amounts to a successive application of the theory for the two- 
part medium. The path, between A and P, consists of three segments whose surface imped- 
ances are Z, Z^ and Zi. The length of the latter two segments are d^ and d^. Using the 




Figure 2. Three-section path on a spherical earth. 

406 



above-mentioned method, it really follows that the expression for the resultant attenuation 
function may be written 



W {d,Z,Z.,Z,)=W(d,Z)-(^-^) -^J^ ______,;, 



ikd\'^Z,-Z pi W{d-a,Z)W(a.Zx) 

/ikdV'^ Z,-Z f .+"2 W{d-a,Z)W'{d,Z,,Z,) , ,._, 



where W{d^Z) is the attenuation function characteristic of propagation from A io B over a 
homogeneous earth of surface impedance Z. W(d—a,Z) and Wia,Zi) are attenuation functions 
for propagation over homogeneous surfaces of surface impedances Z and Zi, respectively. The 
integration variable a can be regarded as a great circle distance measured from B. The 
W^{a,Zi,Z2) occurring in the second integral is the appropriate attenuation function for propa- 
gation over the two-part medium from 5 to a point a on the middle segment. 

For numerical work it is again convenient to introduce the dimensionless parameters used 
for the two-section theory. Thus 

W'{x,q,q2,(li)=W{x,q)+l-.) (qi-q) \. ^' dr. 

\^V Jo [x{x—x)Y^^ 

where 

iq,= (ka/2y^\Z,M, x,= {ka/2y^\d,/a), 

and the other quantities have their usual meaning. From (8) and (9) we may deduce that 

W ix,q,q^^) = W(x,qO+{^-,) (,.-,0j^ [(^-^^T^^ ^''^ 



IM^ ^A_f^Y\ , r^ W{x-x/,Z2)W(x\Z,) , , ,^^. 

= W{x.,q^)M-A {qx-q,) -7^ dx' • (40) 

Vtt^/ Jo \ix—x)xV'^ 



[{x-x')x']'''' 

Of course it would be possible to obtain various series formulas for W^ for a three-section 
path, but in general these would be cumbersome. It is considered to be preferable to work 
directly with (38) and carry out the integration numerically since presumably ^^'(£,^1,^2) 
would be known from the two-section theory. 

There is a rather simple limiting case of (38) which is rather interesting. We take q = q\ 
and X2<^x and Xi. To this approximation, the first integral on the right-hand side of (38) 
vanishes and the integrand in the second integral is essentially a constant over the range Xi 
to Xi+a:;2- Thus 

W{x,q,q,„ci,)^W{x,q)^(^^}j" {q,-q) y^^^^^l ^^-^yn W{x-x,,q)W{x,,q). (41) 

This simple formula has a clear physical interpretation. The second term on the right-hand 
side can be regarded as the field scattered by the strip of width X2. It is proportional to q2—q 
(or Z2—Z) which is the surface impedance contrast and it also is proportional to the attenua- 
tion functions W{x—Xi,q) and W{xi,q) which respectively account for the attenuation from 
A to the strip, and from the strip to jB. Such an equation is analogous to the first Born 
approximation in scattering theory. 

In the foregoing discussion of the three-section theory the antenna heights ha and h^, of 
A and B respectively, may be arbitrary (provided that ha and hb<^a). When the height of 
one terminal is large, such that ^a^l? for example, it is again convenient to define a radiation 
pattern function V\ In this case, 
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ifa 



where 



(g2-g) r^i+"2 F(x-^.g)Ty-(^,gi,g.) ,«; 



where all quantities on the right-hand side have their usual meaning. 

7. Some Extensions of the Theory 

It is of interest to apply the compensation theorem as expressed by (3) to a spherical 
earth in the case when the surface impedances are a function of distance along the great circle 
distance d between A and B. The one-dimensional form of (3) may be written in terms of the 
attenuation functions as follows 

\7rv Jo [x{x—x)Y^^ 

where qi and q are now functions of x (or a) the great circle distance from B. Of course, if 
qi and g are constant, (44) is identical to (7). 

It is believed that (44) could be solved directly by numerical means if q is taken as a con- 
stant corresponding to the unmodified homogeneous earth and qi(x) is some specified function 
of X over the range from to Xi. Equation (44) also permits an approximate treatment of 
transition regions. To illustrate this latter point we consider a three-section path as shown 
in figure 2. Now, however, the middle section is considered to have a surface impedance which 
is a function of a. Thus, ^2 is a function of x over the region from Xi to X1+X2. The attenua- 
tion function for this particular path may be expressed by 

^'"^^ -J-i [xix-x)Y'^ 

(45) 

where W^{x,q,qi) is the appropriate two-section attenuation function and is given explicitly 
by (8). Now on the assumption that X2 is small compared with both Xi and X2 it is permissible 
to replace W^ in the integral in (45) by W{xi,qi). Admittedly, this is a first-order approxima- 
tion but it should be adequate since the relative contribution of the integral is expected to be 
small in any case. To the same approximation W{x—x,q) in the integrand can be replaced 
by W(x—Xi,q). Thus 



where 



T^r-(a-i,g,fe(x),gO^W^X^,g,g.)+(ff''(g2-g) ^^^(^^^^^)]i,, W(a--Xi,g)W(a;i,gO (46) 

q2=— q2{x)dx is the average value of q2{x). 

'^2 Jxi 



8. Appendix 

In the main text of the paper the attenuation function W for a homogeneous earth con- 
sistently appears. It is of interest to see that an integral equation for W in the case ha=hi,=0 
emerges directly from (44). To show this we allow Xi to approach zero in (44) and then we 
imagine the unmodified surface to be a flat, perfectly conducting plane between A and B. 
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Thus, the appropriate vahie W{x,q) is simply exp (ix^ll2). (The quantity x^/12 is the difference 
between the arc and the chord between A and B.) Furthermore, the appropriate value of the 
surface impedance along the arc AB over a conducting plane is easily found to be given by 

or 

Ci{x)^i{x—x)l2. 
Thus 

W^'fe(Zi(a:)) = 6^-^'^^2^('^Y'' r L {x)-i ^^^1 e'''-^'''"'W\x,q,{x))dx. (47) 

Now if gi is a constant throughout the length of the path, W^ in the above equation becomes 
by definition the attenuation function for a homogeneous path. Thus 

This integral equation has been derived by Hufford [1952] who also shows that its solution 
is the residue series 

W{x,ci,)=r-f) ^ f-5 (49) 

where U are roots of 

In a similar manner it may be shown that 

where 

ViX,q^) = -i^^'^ 2 '~T'.. - (51) 

These particular integral equations have also been discussed by Logan [1959]. 
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Selected Abstracts 

Two magneto-ionic phenomena permitting the observation 
of valley minima between the E and F regions in the Arctic, 

J. \\ . Wright, pp. 85-93 of Some ionospheric results obtained 
during the International Geophysical Year; Proc. symp. 
organized by the U RSI I AG 1 committee, Brussels, 1959, edited 
by W. J. G. Beynon (Elsevier Publ. Co., Amsterdam, 1960). 
Two phenomena are described which result from the influence 
of electron collisions on radio wave reflections from the E 
region at high magnetic latitudes. One permits observation 
of the depth of the ''valley" betw^een the E and F regions; 
the other permits the measurement of electron density at a 
fixed and known height in this valley, when it exists. 

The CRPL electron density profile programme: Some features 

and early results, J. W. Wright, pp. 215-220 of Some iono- 
spheric results obtained during the International Geophysical 

Year; Proc. symp. organized by the U US 1 1 AG I committee, 
Brussels, 1959, edited by W. J. G. Beynon {Elsevier Publ. Co., 
Amsterdam, 1960). 

Electron density data and other quantities {hm&xF2, scale 
height, total electron content) provided by reduction of 
vertical soundings to N(h) profiles at CRPL, are described. 
A model of the F region above h,n^xF2 is developed, corre- 
sponding to a Chapman profile of 100 km scale height, and is 
shown to be in agreement with available rocket and moon- 
echo data. Vertical cross sections of the ionosphere across 
the geomagnetic equator are illustrated and discussed. 

Pecularities of the ionosphere in the Far East: Sporadic E 
and F region scatter, E. K. Smith, Jr., and J. W. Finney, 
pp. 182-191 of Some ionospheric results obtained during the 
International Geophysical Year; Proc. sijmp. organized by 

the U RSI I AG I committee, Brussels, 1959, 'edited by W. J. G. 
Beynon {Elsevier Publ. Co., Amsterdam, 1960). 
Through a study of ionosonde data from the world-wide net- 
work of statio.is and of miscellaneous obliciue-incidence 
field-strength measurements made in Japan and the United 
States it was discovered shortly before the beginning of the 
IGY that sporadic E appears considerably more intense in 
the Far East than in similar latitudes in the Western Hemi- 
sphere. Comparisons of sporadic E data are very difficult 
unless made under identical conditions. Therefore an experi- 
ment was designed which consisted of recording transmission 
loss over two matched 50 Mc/s oblique-incidence circuits of 
approximately 800 miles in length, one in the Far East and 
the^ other in the Caribbean. The principal advantage of an 
obhque-incidence circuit over a vertical-incidence one at the 
equivalent frequency (about 11 Mc/s in this case) is that 
although the amount of sporadic E observed would be com- 
parable, the effect of D region absorption in decibels at 
50 Mc/s would be less than at il Mc/s by a factor of about 4.5. 

Studies of scattering phenomena in the equatorial ionosphere 
based upon VHF transmissions across the magnetic equator, 

K. L. Bowles and R. Cohen, pp. 192-195 of Some ionospheric 
results obtained during the International Geophysical Year; 

Proc. symp. organised by the U RSI I AG I committee, Brussels, 
1959, edited by W . J. G. Beynon {Elsevier Publ. Co., Amster- 
dam, 1960). 

Results have been attained pertaining to the phenomena of 
equatorial sporadic E, equatorial spread F, and propagation 
by means of ionospheric radio wave scattering at VHF near 
the magnetic equator. These studies are the results of IGY 
experimentation in South America utilizing 50 Mc/s trans- 
missions crossing the magnetic equator at Huancayo, Peru, 
and over other oblique paths. Additional ionospheric experi- 
ments were also conducted in Huancayo at vertical incidence. 



It has been established that a close relationship exists between 
the magnetic manifestations of the equatorial electrojet cur- 
rent above Huancayo, the occurrence on Huancayo ionograms 
of equatorial sporadic E, and the intensity and fading rate of 
VHF signals propagated by the ionosphere over Huancayo. 
The zone in which E region effects of the equatorial current- 
stream are apparent extends over some 10° of latitude. Also, 
the scatter propagation in the equatorial ionosphere is 
characterized by higher signal-levels and by more rapid 
fading-rates than are observed in the ionosphere at temperate 
latitudes, and signal-strengths nnnain high throughout the 
night. The feasibility of utihzing ionospheric scattering for 
communications in equatorial regions seems c^uite promising. 
Curiously, ionospheric scatter propagation at VHF similar to 
that associated with sporadic E formations in temperate 
latitudes is observed from time to time just to the North and 
South of the magnetic equator, but its occurrence appears to 
be excluded in tlu^ immediate vicinity of that equator. 
Further a remarka])le s(;nsitivity to polarizatioTi was estab- 
lished in HF radar echoes obtained at vertical incidence from 
the diffuse region (bounded by the slant sporadic E and 
extending to the <^-type sporadic E), comprising the typical 
magnetic-ecjuatorial sporadic E configuration appearing on 
Huancayo ionograms. These diffuse echoes were ol)tained 
with an antenna-array polarized magnetic North-South, but 
were virtually (eliminated upon using an East-West polarized 
antenna-array for transmitting and/or receiving. 
Ionospheric propagation by scattering from the F region was 
sought over a transequatorial ]:)ath employing a 2580 km 
transmitter-receiver separation. i*roi)agation by /''-scatter 
was present only about 10 i)ercent of the time, and this was 
during evening hours, wdien enhancements in signal-strength 
above the levels due to Fj region ])ropagation were found to 
occur over this circuit, generally before midnight. A sufficient 
and nearly necessary condition for the occurrence of such 
openings was found to be the ])resence of equatorial spread F 
indications on the Huancayo ionograms. By the variation of 
antenna radiation patterns and by oblique incidence pulse- 
delay measurements, the height of the propagation medium 
sustaining this F-scatter was established, in general, to be 
identifiable with the lowest height of the associated equatorial 
spread F seen on Huancayo ionograms. Thus, the spread-F 
time delays are not the result of an intervening screen below 
the F region that is invisible on the ionograms, as had been 
suggested in certain theoretical explanations. Further, the 
received pulses over the oblique path were but littU* wider 
than those transmitted. This lack of pulse-broadening and 
the narrowing of the equatorial spread F configuration on the 
Huancayo ionogram upon using a high-gain antenna-array 
there at vertical incidence, indicate that the equatorial spread 
F consists of a thin sheet of irregularities. From the pulse 
measurements, the thickness of this sheet is estimated at the 
order of 50 km. The not always perfect comparisons in the 
lowest heights of Huancayo spread F with heights deduced 
from pulse-delay in the F-scatter is interpreted as evidence 
for a limiting geographical extension of a given scattering 
sheet in the F region. Polarization measurements similar to 
those described above prove indirectly that the equatorial 
spread F echoes observed at Huancayo must arrive from the 
magnetic East-West direction; i.e., that orthogonality to the 
magnetic field is imposed as a requirement for obtaining them. 
This is good evidence that the irregularities are elongated in 
the North-South direction, along the earth's magnetic fines 
of force. Thus, the time delays associated with equatorial 
spread F at Huancayo must arise from echoes to the East and 
West of the zenith. From the success in obtaining scattering 
at 50 Mc/s up to heights of 500 km or more, it may be inferred 
that irregularities of scale size 10 metres or smaller, measured 
normal to the magnetic field lines, can exist up to these heights 
at night in the equatorial F region. 
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standard frequencies and time signals from NBS stations 
WWV and WWVH, by Radio Standards Laboratory, NBS 
Misc. Puh. 236 {Dec. 1, 1960). 

Detailed descriptions are given of six technical services 
broadcast by National Bureau of Standards radio stations 
WWV and WWVH. The services include 1, standard radio 
frequencies; 2, standard audio frequencies; 3, standard time 
intervals; 4, standard musical pitch; 5, time signals; and 6, 
radio propagation forecasts. Other domestic and foreign 
standard frequency and time signal broadcasts are tabulated. 

Near infrared atmospheric transmission to solar radiation, 

D. M. Gates, /. Opt. Soc. Am, 50, No. 12, 1299-1304 {Dec. 
1960) . 

Near infrared solar spectrum observations taken on October 
15, 1954 with a double-pass NaCl prism spectrometer have 
been analyzed for transmission coefficients for the ''selective" 
absorption factor and for the ''continuum" factor. The 
analysis was carried out for 59 wavelength positions between 
0.872 and 2.537m. The monochromatic data lit well the law 
lnT=ci{w)^^^ where w is the amount of water vapor in the 
optical path. The coefficient ci is given as a continuous 
function of the wavelength. A coefficient of extinction for 
the "continuum" factor is also given. 

Limitations of radiosonde punch-card records for radio- 
meteorological studies, B. R. Bean and B. A. Cahoon, 
/. Geophys. Research 66, No. 1, 328-331 {Jan. 1961). 
Instrumental and data reduction restrictions affect the con- 
version of radiosonde data into radio refractive index profiles. 
It is found that over half of the elevated radio ducts reported 
by the observer will be missed due solely to the standard 
criteria established for rejection of data points in the process 
of preparing punch-card records. 

An extensive collection of radiosonde data containing all 
observed levels is now available from some seventy locations 
on and near the North American Continent. 

VLF phase perturbation associated with meteor shower 
ionization, C. J. Chilton, /. Geophys. Research 66, No. 2, 
379-383 {Feb. 1961). 

Anomalies associated with the Lyrid, 5-Aquarid, and Perseid 
meteor showers were observed on the daily phase variation 
of the 16 kc/s transmission from Rugby, England, to Boulder, 
Colorado. The effective reflection heights estimated from 
these anomalies are 80, 81, and 83 km respectively. 

Microwave spectroscopy— Atomic frequency standards, J. M. 

Richardson, Encyclopedia of Spectroscopy, edited hy G. L. 
Clark {Reinhold Publ. Corp., New York, N.Y., 1960). 
The topic of atomic frequency standards is surveyed in a 
general way with emphasis on principles and current results. 
Atomic transitions satisfy very well certain requirements of 
a good standard; namely, constancy, precision, renewability, 
and convenience. Technique of observation, results, and 
limitations of several types of atomic frequency standards 
are described. These are the gas absorption, atomic beam, 
maser, and optical pumping methods. The atomic resonance 
may be used as a resonator or as a frequency controlling 
element in a control mechanism. Useful applications are 
improvement in the time scale, analysis of power spectra and 
line shapes, and tests of certain postulates of special and 
general relativity. 

A theoretical study of sporadic-E structure in the light of 
radio measurements, K. Tao, NBS TN87 {PB161588) 
{$1.25). 

The theoretical aspects of the mechanisms of sporadic-jB/ 
reflections are described from both the standpoint of a thin 
layer and a scattering model. For the thin layer model, thin 
dielectric layers which have various distributions of electron 
density are considered. It is also pointed out that the scatter- 
ing theory for which an autocorrelation function of the 
fluctuation of electron density is given by modified Bessel 
functions of the fourth through the seventh order is an 
available model for sporadic-J^ scatter. Moreover blobs of 
ionization which have a horizontal scale of the order of 200 m 
and a vertical scale of about 50 m are considered for sporadic-^ 
scatter. The frequency and distance dependences of the 
oblique VHF propagation by means of the sporadic-^E/ layer 
are discussed by comparing the theoretical results with 
experimental evidence. 



The NBS meteor-burst propagation project — a progress report, 

C. E. Hornback, L. D. Brevfogle, and G. R. Sugar, NBS 
TN86 {PB 16 1587) {I960) $L25. 

This report briefly describes a meteor-burst propagation- 
study program at NBS-Boulder Laboratories and presents 
some of the preliminary analysis results. Observations have 
been made with scaled systems over three different paths 
(Long Branch-Table Mesa, Norman-Fargo, and Barrow- 
Kenai) at frequencies of 30, 50, and 74 Mc/s. The recorded 
data is processed by a combination of manual and automatic 
methods. The preliminary results show about a 10-db diurnal 
variation in threshold for a constant duty-cycle. Thresholds 
for a constant duty-cycle were observed to have an approxi- 
mate frequency dependence relative to 30 Mc/s of 15 db 
lower for 50 Mc/s and 30 db lower for 74 Mc/s. There was 
no statistically-significant difference observed in the 
occurrence of meteor-bursts from a Poisson distribution. 

Oblique incidence receiving antenna array for a relative 
ionospheric opacity meter, A. C. Wilson, NBS TN78 

{PB161579) {Nov. 1960) 50 cents. 

Experimental measurements incidental to the design of an 
antenna for a relative ionospheric opacity meter (RIO 
Meter) are described. 

The frequency of operation is 50 Mc. The antenna require- 
ments are that the main lobe of the antenna is directed at 
23° above the horizon, the half-power beam-width in the 
vertical plane does not exceed 10°, the minimum front-to-back 
ratio is at least 13 decibels, and the side-lobe levels are at 
least 10 decibels below the maximum response in bo*h the 
E- and H-planes. Since the antenna is for use in Alaska, it 
is to be of simple design and physically able to withstand any 
anticipated wind and ice loads. 

The final antenna design is an array of three stacked horizontal 
dipoles with two optimally spaced reflectors behind each 
dipole to obtain the required directivity. The narrow main 
lobe of the antenna directed at an angle of 23° above the 
horizon is obtained by properly spaced and phased dipoles 
above the ground. The half-power beamwidth in the vertical 
plane is computed to be 7)i°. The front-to-back ratio over 
the rear 180° sector is not less than 20 decibels, and the 
half -power beamwidth in the E-plane is 74°. 
Two complete receiving antenna arrays were constructed, 
adjusted, and tested. These antennas were installed in 
Alaska where they are now in use by an auroral transmission 
loss project. 

A study of auroral coruscations, W. H. Campbell and M. H. 
Rees, J. Geophys, Research 66, No. 1, 41-55 {Jan. 1961). 
Short period variations in the N2+(0, 0) auroral emission 
band represent 5 per cent of the total light in the 3914A 
region. These quasi-periodic coruscations have a dominant 
period of 6 to 10 seconds. They attain a maximum amplitude 
in the predawn hours and are closely related to magnetic 
field micropulsations and ionospheric absorption of cosmic 
noise. Spectroscopic triangulation showed that the variations 
originate in the E region of the ionosphere. The profile of 
the electron density associated with the aurora was found to 
have a maximum value of 1.1X1<0^ electrons per cubic 
centimeter at 98 kilometers. 

Reports on CSAGI Disciplines, Part Ai, The Airglow, hy 
F. E. Roach. Annals of the International Geophysical Year, 

X 134-137 {Pergamon Press, New York,^ 1960). 

Hourly zenith intensities of the night airglow in rayleighs for 

IGY and IGC. Results from 28 stations. 

Propagation of electromagnetic pulses in a homogeneous 
conducting earth, J. R. Wait. Appl.Sci. Research, Section B, 
8, 213-253 {I960). 

A general analysis for the electromagnetic response of con- 
ducting media due to pulse excitation is presented. The 
treatment is based on the Laplace transform theory. First, 
a survey of the field is made and the limitations and scope of 
the previous work are pointed out. The theory of propaga- 
tion of a plane wave pulse in a conducting and homogeneous 
medium of infinite extent is then reviewed. The form of these 
results enable one to evaluate the relative importance of the 
conductivity and the dielectric constant. It is indicated, for 
sufficiently large times in the transient response, that displace- 
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meiit currents may be safely neglected for sea water and for 
most geological media. Under this assumption, the wave- 
form of the electric field in a conducting medium is illustrated 
for the case where the source is an electric dipole energized by 
a step-function current. Results are also presented for 
exponential and bell-shaped source functions. The pulse 
shape of the field components is profoundly modified as they 
propagate through the medium. It is suggested that this 
proi)erty may be utilized in measuring distances in the earth's 
crust. The more difficult problem of propagation in non- 
infinite conducting media is also considered. To account for 
the presence of the interface in a conducting half space (i.e., 
homogeneous flat ground), a rather involved analytical ex- 
pression for the transient fields is required. Certain special 
cases, such as a horizontal electric dipole at the interface, are 
illustrated by numerical results. The transient excitation of 
a wire loop lying on the surface of a homogeneous ground is 
also considered. Finally, transient coupling between pairs of 
parallel insulated wires grounded at their end points is treated 
as an extension of the earlier results. 

Ionospheric absorption at times of auroral and magnetic 
pulsations, W. H. Campbell and H. Leinbach, /. Geophys. 
Research 66, No. 1, 25-35 {Jan. 1961). 

A study in March and April 1960, showed variations in the 
auroral zone ionospheric absorption of cosmic noise to be 
closely related to magnetic field micropulsations and short 
period coruscations of X 3914. At times of polar-cap type 
absorption, magnetic micropulsation amplitudes were dimin- 
ished. Auroral ionization in the E region, estimated from a 
particular luminosity-height profile, accounted for 50 per 
cent, at least, of the cosmic noise absorption. 

Tests for regression coefficients when errors are correlated, 

M. M. Siddiqui. Annals of Math. Stat. 31, No. ^ 929-938 
{Dec. 1960). 

In a previous paper the covariances of least-squares estimates 
of regression coefficients and the expected value of the 
estimate of residual variance were investigated when the 
errors are assumed to be correlated. In this paper we will 
investigate the distribution of the usual test statistics for 
regression coefficients under the same assumptions. Applica- 
tions of the theory to the cases of testing a single sample mean, 
the difference between the means of two samples, the coeffi- 
cients in a linear trend and in regression on trigonometric func- 
tions will be discussed in some detail under an assumed covari- 
ance matrix for errors. 

Climatic charts and data of the radio refractive index for the 
United States and the World, B. II. Bean, J. D. Horn, and 
A. M. Ozanich, Jr., N BS Mono. 22 {Nov. 25, 1960) $2.00. 
The radio refractive index of air, n=l-\-NX 10~^, is a function 
of atmospheric pressure, temperature, and humidity and 
varies in a systematic fashion with climate. 
Included in this Monograph is a compilation of refractive 
index data. Data listings made up of observations from 45 
U.S. surface weather stations for 2-hour intervals over an 
8-year period are given. Mean values, maxima, minima, and 
standard deviations of the refractive index have been cal- 
culated and tabulated for these observations. Additionally, 
mean vertical profiles of the refractive index have been pre- 
pared for 43 U.S. upper air sounding stations from long-term 
means of pressure, temperature, and humidity. 
Earlier studies of refractive index climate are assimilated and 
put into perspective. One such study is an extensive analysis 
and mapping of the refractive index climate of the United 
States. A worldwide radio refractive index climatology is 
developed based upon monthly mean observations of pressure, 
temperature, and humidity. 

An important finding of these climatological investigations is 
the strong correlation of N with height. A reduced-to-sea- 
level value of the index, termed No, is used to eliminate this 
systematic height dependence. The surface value of N, Ns, 
may be estimated four to five times more accurately from 
charts of A^o than from similar-sized charts of Na itself. 
From climatic charts of No, Ns niay be estimated at any given 
location in the United States throughout the day during every 
season. In addition detailed annual and diurnal cycles, as 
well as 8-year cumulative probability distributions, are given 
for 12 representative U.S. stations. 



On a worldwide basis, charts of mean No are presented for 
both summer and winter season. 

Diffractive corrections to the geometrical optics of low- 
frequency propagation, J. R. Wait. Electromagnetic Wave 
Propagation (International Conference Sponsored by the 
Postal and Telecommunications Group of the Brussels Uni- 
versal Exhibition), edited by M. Desirant and J. L. Michiels, 
pp. 87-101 {Academic Press, New York, N.Y., 1960). 
The influence of caustics in low frequency radio propagation 
is studied. Correction factors to be applied to the geometrical 
optical representations are presented in the form of graphs. 
Particular attention is paid to the treatment of the caustic 
w^hich occurs at the geometrical horizon of the first hop sky 
wave. It is shown that the apparent infinity in the con- 
vergence coefficient can be removed when the wave nature of 
the problem is considered. Antipodal phenomena are also 
discussed. 

Note. In this paper the publisher deleted the captions 
for Figures 7a, 7b, 8a, and 8b. These should read: 

Fig. 7a — The factor \F\ for sea water, e = 80eo, a = 5 mhos/ 
meter, a = 4/3X4360 km. 

Fig. 7b — The phase, of the factor F, for sea water. 

Fig. 8a — The factor |/''| for average land, e= 15€o, <r = 5X 10~^ 
mhos/meter, a = 4/3X4360 km. 

Fig. 8b — The phase A, of the factor F, for average land. 

Furthermore, in these four figures the abscissae should be 
lal^elled "Angle below horizon (0 — 0c)/2." Also, in equation 
(25) the factor (/'') is defined by 

where |F| and A are the quantities plotted in Figures 7 and 8. 

Propagation of electromagnetic waves along a thin plasma 
sheet, ,1. U. Wait, Can. J. PInjs. 38, 15Sh'~15!L^, {Dec 1960). 
It is shown that a thin ionized sheet will support a trapped 
surface wave. The effect of a constant and uniform magnetic 
field is to modify the phase velocity and polarization of the 
surface wave. The essential features are illustrated by nu- 
merical results for selected values of the electron density, 
collision frequency, and gyro freqiuuicy. The effect of locat- 
ing the plasma sheet near and parallel to a conducting plane 
is also considered. In this situation other modes of a wave- 
guide type are possible in addition to the surface wave. 

Ionospheric mapping by numerical methods, W. li. Jones and 
H. M. (lallet, Tekcomm. J., Journal LIT, No. 12, 260e-264e 
{Dec. 1960). 

One of the most difficult problems encountered in forecasting 
conditions for long-distance radio propagation is the mapping 
of ionospheric characteristics. A satisfactory method has 
been developed for mapping the continuous time variations on 
a world-wide basis of any ionosplieric quantity entirely by 
numerical methods and electronic computing. 
As an illustration, the method is applied to mapping the 
monthly median values of the critical frequency of the F2 
layer (foF2) for December 1957. Input to the computer con- 
sists of the data as they are tabulated at 109 available sta- 
tions. An analysis is performed l)y progressively fitting series 
of orthogonal functions in three coordinates: time, latitude, 
and longitude. Two main difficulties have been overcome: 
(1) the data were irregularly positioned in the two space 
dimensions, and (2) they were affected by statistical noise. 
The end product of the analysis is a numerical map, or mathe- 
matical function, defined by a table of coefficients, small in 
comparison with the number of input data. 
Many applications can be made using the table of coefficients. 
For example, a variety of graphical representations and sta- 
tistical tests for accuracy may be produced, as well as fore- 
casting ionospheric conditions using a suitable solar index. 
The numerical methods are general enough to be applied to 
any geophysical quantity much as meteorological or geomag- 
netic characteristics varying continuously with time on a 
world-wide basis. 
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other NBS Publications 

Journal of Research, Vol. 65A, No. 3, May-June 1961. 
70 cents. 

International practical temperature scale of 1948. Text 

revision of 1960, H. F. Stimson. 
Evaluation of the NBS unit of resistance based on a com- 
putable capacitor, R. D. Cutkosky. 
Wavelengths and intensities in the first spectrum of bromine, 

2000 to 13000 A, J. L. Tech and C. H. Corliss. 
Torsional resonance vibrations of uniform bars of square 

cross section, W. E. Tefft and S. Spinner. 
Infrared studies of aragonite, calcite, and vaterite type 

structures in the borates, carbonates, and nitrates, C. E. 

Weir and E. R. Lippincott. 
Dielectric properties of polyamides: polyhexamethylene 

adipamide and polyhexamethylene sebacamide, A. J. 

Curtis. 
Heat of formation of calcium aluminate monocarbonate at 

25 °C, H. A. Berman and E. S. Newman. 
Thermodynamic constants for association of isomeric chloro- 

benzoic and toluic acids with 1,3-diphenylguanidine in 

benzene, M. M. Davis and H. B. Hetzer. 
Heats of combustion and formation of trimethylborane, 

triethylborane, and tri-n-butylborane, W. H. Johnson, 

M. V. Kilday, and E. J. Prosen. 
Pyrolvsis of linear copolymers of ethylene and propylene, 

S. Straus and L. A. Wall. 
Pyrolysis of fluorocarbon polymers, L. A. Wall and S. Straus. 
Preparation of fluoro- and bromofluoroaryl compounds by 

copyrolysis of bromofluoroalkanes. L. A. Wall, J. E. 

Fearn, W. J. Pummer, and R. E. Lowry. 
Thermal stability of polydivinylbenzene and of copolymers of 

styrene with divinylbenzene and with tribinylbenzene, 

S. Straus and S. L. Madorsky. 
Conformations of the pyranoid sugars. IV. Infrared absorp- 
tion spectra of some fully acetylated pyranoses, R. S. Tip- 
son and H. S. Isbell. 
A standard for the measurement of the pT{ of blood and other 

physiological media, V. E. Bower, M. Paabo, and R. G. 

Bates. 

Journal of Research, Vol. 65B, No. 2, April-June 1961. 
75 cents. 

Optimal approximation for functions prescribed at equally 

spaced points, H. F. Weinberger. 
Truncations in the method of intermediate problems for 

lower bounds to eigenvalues, N. W. Bazley and D. W. Fox. 
Comparison theorems for symmetric functions of character- 
istic roots, M. Marcus. 
Some properties of the empirical distribution function of a 

random process, M. M. Siddiqui. 
Another extension of Heinz's inequality, M. Marcus. 
Mean motions in conditionally periodic separable systems, 

J. P. Vinti. 
Some boundary value problems involving plasma media, 

J. R. Wait. 
A new decomposition formula in the theory of elasticity, 

J. H. Bramble and L. E. Payne. 
Pointwise bounds in the Cauchy problem of elastic plates, 

L. E. Payne. 

Journal of Research, Vol. 65 C, No. 2, April- June 1931. 75 
cents. 

An experimental study concerning the pressurization and 

stratification of liquid hydrogen, A. F. Schmidt, J. R. 

Purcell, W. A. Wilson, and R. V. Smith. 
Temperature dependence of elastic constants of some 

cermet specimens, S. Spinner. 
Analog simulation of zone melting, H. L. Mason. 
Residual losses in a guard-ring micrometer-electrode holder 

for solid-disk dielectric specimens, A. H. Scott and W. P. 

Harris. 
A bolometer mount efficiency measurement technique, G. F. 

Engen. 



Telescope for measurement of optic angle of mica, S. Ruthberg. 

An automatic fringe counting interferometer for use in the 

calibration of line scales, H. D. Cook and L. A. Marzetta. 



Units of weight and measure (United States customary and 

metric) definitions and tables of equivalents, NBS Misc. 

Publ. 233 (1960) 40 cents. 
Report of the 45th national conference on weights and 

measures 1960, NBS Misc. Publ. 235 (1960) 75 cents. 
Measurements of absorbed dose of neutrons, and of mixtures 

of neutrons and gamma rays, NBS Handb. H75 (1961) 

35 cents. 
Medical X-ray protection up to three milhon volts, NBS 

Handb. H76 (1961) 25 cents. 
Precision measurement and calibration, S. F. Booth: Elec- 
tricity and electronics, NBS Handb. 77, Vol. I (1961) 

$6.00; Heat and mechanics, NBS Handb. 77, Vol. II (1961) 

$6.75; Optics, metrologv, and radiation, NBS Handb. 77, 

Vol. Ill (1961) $7.00. 
Report of the International Commission on radiological units 

and measurements (ICRU) 1959, NBS Handb. H78 (1961) 

65 cents. 
Quarterly radio noise data — June, July, August 1960, W. Q. 

Crichlow, R. D. Disnev, and M. A. Jenkins, NBS TN18-7 

(PB151377-7) (1960) $1.75. 
VHF radio propagation data for Cedar Rapids-Sterling, 

Anchorage-Barrow, and Fargo-Churchill test paths April 

1951 through June 1958, G. R. Sugar and K. W. SulHvan. 

NBS TN79 (PB161580) (1960) $4.00. 
Bibliography of tropospheric radio wave scattering, R. L. 

Abbott, NBS TN60 (PB161581) (1960) $2.25. 
A survey of spread-/, F. N. Glover, NBS TN82 (PB161583) 

(1960) $1.75. 
Some ceramic dielectrics with a very low temperature coeffi- 
cient of capacitance, S. MarzuUo and E. N. Bunting, 

J. Am. Ceram. Soc. 43, No. 11, 609, (Nov. 1960). 
Vibration-rotation bands of N2O, E. D Tidwell and E. K. 

Plyler, J. Opt. Soc. Am. 50, No. 12, 717-720 (Dec. 1960). 
Nuclear magnetic resonance in tantalum metal, J. I. Budnick 

and L. H. Bennett, J. Phys. Chem. Sohds 16, No. 1/2, 37-38 

(Jan. 12, 1960). 
The sample, the procedure, and the laboratory, W. J. Youden, 

Anal. Chem. ^2, No. 13, 23A-37A (Dec. 1960). 
S-quinolinol precipitation of the elements, J. I. Hoffman, 

Chem. Anal. 49, No. 4, 126 (Dec. 1960). 
Dipole moments of hvdrocarbons, D. R. Lide, Jr., J. Chem. 

Phys. 33, No. 6, 1879 (Dec. 1960). 
Collection of ions produced by alpha particles in air, Z. Bay 

and H. H. Seliger, Phys. Rev. 120, No. 1, 141-143 (Oct. 1, 

1960). 
A test of a procedure for easy estimation of representative 

monthly electron density profiles for the ionosphere, 

J. W. Wright, J. Geophys. Research 65, No. 10, 3215-3217 

(Oct. 1960). 
A method for efficiently providing low temperature liquids on 

a large scale to an accelerated experimental program, 

A. K. Stober, Proc. 10th Intern. Congress of Refrigeration, 

Copenhagen, Denmark 1, No. 1, 17-18 (Pergamon Press 

Ltd., London, England, 1959). 
Seasonal variations in the twilight enhancement of [01] 5577, 

L. R. Megill, P. M. Jamnick and J. E. Cruz, J. Atmospheric 

and Terrest. Phys. 18, 309-314 (Aug. 1960). 
Synthetic mica, A. V. Valkenburg, Book, Encyclopedia of 

Chemical Techniques, p. 480-487 (Interscience Encyclo- 
pedia, Inc., New York, N.Y., 1960). 
Nonresonant microwave absorption and relaxation frequency 

at elevated pressures, A. A. Marvott and G. Birnbaum, 

J. Phys. Chem. 64, No. 11, 1778-1780 (Nov. 1960). 
Precipitation of the elements on addition of aqueous ammonia 

to their clear solution, J. I. Hoffman, Chem. -Anal. 49, 

No. 3, 94 (Sept. I960). 
Low scatter high current gas target for D-D neutrons, A. C. B. 

Richardson, Rev. Sci. Inst. 31, No. 11, 1202-1203 (Nov. 

1960). 
The heat of combustion of dicyanoacetylene, G. T. Armstrong 

and S. Marantz, J. Phys. Chem. 64, 1776-1777 (1960). 
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Evaluation of resistance strain gages at elevated temperatures, 

R. L. Bloss, Am. Soc. Testing Materials Proc. 1, No. 1, 

9-15 (Jan. 1961). 
Errors in dielectric measurements due to a sample insertioxi 

hole in a cavitv, A. J. Estine, and H. E. Bussev, IRE Trans. 

Microwave Theory and Tech. MTT-8, No. 6, 650-653 

(Nov. 1960). 
Ensemble method in the theory of irreversibility, R. Zwanzig, 

J. Chem. Phys. 33, No. 5, 1338-1341 (Nov. i960). 
Determination of pentosans. Interlaboratory comparison of 

the aniline acetate, orcinol, and bromination methods, W. K. 

Wilson and J. Mandel, Tappi 43, No. 12, 998-1004 (Dec. 

1960). 
A research for geomagnetic singular days, C. Warwick, T. 

Pohrte, and N. MacDonald, J. Geophys. Research 65, No. 

9, 3013-3015 (Sept. 1960). 
Letter to editor of Revue Des Materiaux (first), E. S. Newman 

and H. A. Berman, Revue Des Materiaux C, No. 540, 231 

(Dec. 1, 1960). 
Geomagnetic storms and the space around the earth, S. Chap- 
man, Nature 187, No. 4740, 824-827 (Sept. 3, 1960). 
On the theory of the critical point of a simple fluid, M. S. 

Green, J. Chem. Phys. 33, No. 5, 1403-1409 (Nov. 1960). 
Some implications of slant-/i/s, E. K. Smith and R. W. Knecht, 

Polar Atmosphere Symp., Part 2, Ionospheric Section, Oslo, 

Norway, Julv 2-8, 1956, p. 195-204 (Pergamon Press, Inc., 

New York, N.Y., 1957). 
Mismatch errors in microwave phase shift measurements, 

G. E. Schafer, IRE Trans. Microwave Theory and Tech. 
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